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Abstract. We compute the waves propagating on a compact 3-manifold of constant positive 
curvature with a non trivial topology: the Poincare dodecahedral space that is a plausible model 
of multi-connected universe. We transform the Cauchy problem in a mixed problem posed on a 
fundamental domain determined by the quaternionic calculus. We adopt a variational approach 
using a space of finite elements that is invariant under the action of the binary icosahedral group. 
The computation of the transient waves is validated with their spectral analysis by computing a lot 
I of eigenvalues of the Laplace-Beltrami operator. 
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I- Introduction. 

Some fundamental open questions regarding the nature of the universe concern its geometry and 
topology. They are the subject of many articles ([5], [l2], [13], [IT], [20] for example) in which 
models are compared with observations which are more accurate over the years. About geometry, 
'q_^ after nine years of investigation, the data by WMAP (Wilkinson Microwave Anisotropy Probe) 
mission have provided strong evidence suggesting that the universe is nearly flat, with the ratio of 
its total matter-energy density to the critical value very close to one [231, but without fixing the sign 
j3 of its curvature. The Poincare dodecahedral space (PDS) is a plausible multi-connected universe 
model with positive spatial curvature. More precisely, PDS is the quotient of the unit 3-sphere 
5^ by the group I* of covering isometries. The local geometry of the universe is described by a 
T-H Friedmann-Lemaitre metric, so PDS (denoted by K in equations) is endowed with the spherical 
metric ds^ induced by ds^a . It is a three dimensional C°° compact manifold, without boundary, 
^y-^ Such a manifold can conveniently be converted to a description as a fundamental domain. It will 

have the shape of a dodecahedron, with pairs of faces identified. 
^ Many authors have studied PDS ( [I],[3], [7], [H], [14j, [15j, [H], [21]). A stationnary approach 
CN has been used in all these articles: the exact knowledge of the eigenvalues of the Laplace-Beltrami 
operator Ak on PDS allows them to predict some results (for example the cosmic microwave 
I background temperature anisotropics), and to compare with the WMAP observations. 

In this paper we adopt the time approach previously developed in |2j, where we studied a toy 
model of hyperbolic universe. Among several advantages, this method could be extended to in- 
^ vestigate numerically the non-linear dynamics. The construction of PDS is detailled in the next 
section. In the third part we get a practical description of the particular fundamental domain J- 
that contains (1, 0, 0, 0) of M^, and also we get a numerical description of the projection of T in 
M^. We mainly use the quaternionic calculus. Next we consider the scalar wave operator — Ak 
on the dodecahedral universe x K. Then, we compute the solutions of the wave equation in the 
time domain, by using a variational method and a discretization with finite elements. The domain 
of calculus is J-^, therefore the initial Cauchy problem on the manifold without boundary K, be- 
comes a mixed problem on J^y with suitable periodic boundary conditions on dJ-'v These boundary 
constraints are implemented in the choice of the basis of finite elements. We validate our results 
by performing a Fourier analysis of the transient waves that allows to find a lot of eigenvalues of 
the Laplace-Beltrami operator Ak on PDS. 
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II. The Poincare Dodecahedral Space. 

To be able to perform computations on PDS we need to describe it accurately. This section 
recalls the properties of S^, I* and that we need to know. The 3-sphere is the submanifold 

of four-dimensional Euclidean space such that Xq + + + X3 = 1 . The comoving spatial distance 
d between any two points x and y on 5^ is given by: 

d{x,y) = arccos[x*yj] 

We also use the parametrisation: 

Xq = cosx 
xi = sin X sin sin 99 
X2 = sin X sin cos (/? 
X3 = sin X cos 

with 0<x^'''", O<0<7r and < < 27r. These parametrisation leads to a convenient way to 
visualize S^: two balls in glued together by their boundary. The 9 and 93 coordinates are the 
standard ones for 5^. For the first ball x runs from at the center through | at the surface. For 
the second ball x runs from ^ at the surface through n at the center. We use a projection p of 
in M^: 

p: ^ M3 

(xo,Xi,X2,X3) 1-^ (xi,X2,X3) 

That is to say, to represent a point of coordinates (a;i){j=o,i,...,3} of '^^j consider only the coordi- 
nates (2;j){j=i^...^3}. and discard xq. The two points of coordinates (xq, xi, X2, X3) and (— xq, xi, X2, X3) 
have the same coordinates (xi,X2,X3); they are represented by the same point. In the sequel, for 
any ^4 C 5^, we use Ay to denote the visualization in of A, that is A^ := p{A). 

To define a discrete fixed-point free subgroup F C SO (4) of isometrics of S^, one makes use of 
the fact that the unit 3-sphere is identified with the multiplicative group Hi of unit quaternions 
by 

X = (xo, xi, X2, X3) <S=^> q = xqI + xii + X23 + xsk. 

The four basic quaternions {1, i, j, k} of the set EI of all quaternions, satisfy the multiplication rules 
= = k^ = — 1, and k = ij = — ji. They commute with every real number. The norm of q 
is defined by \\q\\ := Xq -|- xf -|- x^ + Xg. The conjugate g of g = al -|- 6i -|- cj -|- dk is defined by 
g := al — 6i — cj — dk. 

In this paper, we are interested in the subgroup I* of SO (4) called the binary icosahedral group 
([19], [22], [16], [7]). It is a two sheeted covering of the icosahedral group Z C SO(3) consisting of all 
orientation-preserving symmetries of a regular icosahedron. Indeed, for any q in the multiplicative 
group of unit quaternions Hi we consider pg : H — )• H defined by Pq{q') := qq'q~^- It fixes the 
identity quaternion 1, so in effect its action is confined to the equatorial 2-sphere spanned by 
the remaining basis quaternions (i,j,k). Thus, by viewing as the space of pure imaginary 
quaternions, subspace of H with basis {i,j,k}, we get a rotation in . So Pg belongs to SO (3). 
The map vr : Hi — )• S0(3) defined by Tr{q){q') = qq'q~^ is a homomorphism of H onto SO(3). vr is 
a two to one homomorphism as 7r(g) = 7r(— g). By definition we have I* := 7r^^(X). The order of 
X is 60, so the order of X* = tt~^(I) is 120. X* contains only right-handed Clifford translations 7, 
in other words 7 acts on an arbitrary unit quaternion g G 5^ by left multiplication and translates 
all points qi,q2 £ by the same distance Xi d{qi,^qi) = ^((72,7^2) = X- The right-handed 
Clifford translations act as right-handed cork screw fixed-point free rotations of . 
Let take 7 = al -|- 6i -|- cj -|- dk an element of X* (with arccos a = x)i and q = xqI + xii -|- X2j + X3k 
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in 5^. We can express 7 as a matrix and we have: 



7 q 



a —b — c 

b a —d 

c d a 

d —c b 




X* is generated by two isometries s and 7. Denoting by o" = (1 + ^/b)/2 the golden number we 
have: 

X* = (5,7|(57)' = s' = 7'> 

= ibi, ibj, ibk, -(±1 lb i lb j lb k), -(01 ib i ib — j ib dk) {with even permutations) 

We see that every element of I* is a right-handed Clifford translation with x equal to 0, vr, ^, |, 
if' f ' T"' T"' T"" deduce that ^ is the smallest non-zero translation distance. For example 
X* is generated bys = ^(1-bi-bj-bk) and 7 = |1 -b i^j ~ 5 k. s is a right-handed Clifford 
translation with x equal to |. Also, ps is a rotation in M^, its angle is ^ and its axis is directed 
by the vector Ij !)• 7 is a Clifford translation with x equal to |. Moreover p^ is a rotation 

in M^, its angle is ^ and its axis is directed by the vector ^^g!— (0, i, — 1). 

We now consider S^/I* the quotient of the 3-sphere under the action of the discrete fixed- 
point free subgroup I* of isometries of S^, with I* acting by left multiplication. This is the 
Poincare dodecahedral space (PDS) ([E], [22], [16], [7]). To perform the computations of the waves 
on it is very useful to represent it by a fundamental domain C and an equivalence 

relation ~ such that 

(II.l) = TI ^ . 

J- is such that: 

• 53 = Utex.t(-F), 

• Vt G X*, Vt' G X*, t(T) n t'(J')= 0. 

is a regular spherical dodecahedron (dual of a regular icosahedron), and ~ is obtained by identi- 
fying any pentagonal face of with its opposite face, after rotating by ^ in the clockwise direction 
around the outgoing axis orthogonal to this last face. 120 such spherical dodecahedra tile the 
3-sphere in the pattern of a regular 120-cell (see figure [s]). 

III. Fundamental domain T of PDS, and its visualization T^. 

We have choosen the unique fundamental domain that contains (1, 0, 0, 0) to simplify calculations 
and visualization. In the following it will be denoted T . In order to perform computations we need 
to know the coordinates of all points in T ^ and the equations of its edges and faces. We also deduce 
the characteristics of J-^ := p{J~), which is our domain of calculus and visualization. 

Proposition III.l. The unique fundamental domain that contains (1, 0, 0, 0) is the geodesic convex 
hullmS^ of these 20 vertices: ^(a^, 0, ±1), ^(a^, 0, ±1, ±^), ^(^r^, ±i, ±i, ±i) and 

^(a^±l,±^,0). 

Anyone of the 12 faces Fi of F is a regular pentagon in included in 

{(xo,x,y,2;) G 5^, aix + hy + CiZ = ^} , 

with {ai,bi,Ci) equal to (±^,±1,0), up to an even permutation. 

The set of all barycenters in of the vertices of Fi, denoted by F^, is included in a 2-plan of 

2 

equations xq = and aix + biy + CiZ = 
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Proof. J- is the geodesic convex hull in of its 20 vertices, so we begin to search the coordinates 
of the vertices. As the shortest translation distance of the elements of I* is ^, the shortest trans- 
lation distance in PDS is ^. Two opposite faces must be the image from each other by a Clifford 
translation with x equal to |. I* has 12 such elements denoted by gi in the following parts: 

-(ctI + i lb — j lb k) (with even permutations of the tree last coordinates). 
2 a 

So the isometrics of SO{3) denoted by pg. leave invariant p{J^) := J^v Furthermore (0, 0, 0) belongs 
to because (1,0,0,0) belongs to J". So OBi := ^^(0, ±1, ±1), with even permutations, are 
the orthogonal axes to each pair of opposite faces of J-"^. Bi are the vertices of an icosahedron 
that is the dual of a regular dodecahedron whose vertices are the barycenters of three equidistant 
vertices of the icosahedron. We get a pentagonal face having OBi as symetry axis by finding the 
five vertices that are at the same minimal distance from B^. Among these five points, two vertices 
are adjacent if they are at a distance inferior or egal from the other vertices of this face. We deduce 
that, in M^, the last three coordinates of the vertices of J- could be: 

(ib-cr, zb-fj, ib-o"), or (0,±-a'^ ,±-) (with even permutations). 
6 6 6 6 6 

So the coordinates of the 20 vertices Cj belonging to and being the vertices of a regular dodec- 
ahedron that has the same symetry axes are of the form: 

v /l-3(A§ )'l +X{±lai±lai±lak), 
or Y^l — -b A(Oi =b g(T^j lb gk) {with even permutations), 

with A G M such that two opposite faces can be the image from each other by a Clifford translation 
with X equal to |. Consider a face F oi F such that ^i — j)isa symetry axis of J-^ orthogonal 
to the induced face of Fy. Its adjacent vertices satisfy: 

Ci = Vl-3(A|)'l +A(i(-cTi-aj + ak)), C2 = Vl-3(Af)'l + A(| (-a^i - j)), 
Cs = Vl-3(A|)'l + A(i (-ai - - ak)), C4 = Vl-3(Af)'l + A(i (-a^j - k)), 
C5 = Vl-3(A|)'l + K\ (-cT^j + k)). 
And the opposite face having the same symetry axis has the following adjacent vertices: 

Q = v/l-3(Af)'l + \{\ (ai + - ak)), C7 = 7 l-3(A§)^ l + A(| (a^i + j)), 
= Vl-3(Af)'l + \{\ (ai + + ak)), C9 = yi - 3 (AD'l + X{\ (a^j + k)), 
Cio = Vl-3(A|)'l + \{\ - k)). 

We are searching A such that Ci has Cg as image by the Clifford translation ^(cl + ^i+j). So 
the spherical distance between Ci and Cs is equal to |. We get A^ = |(5 — 3it) = Due to the 
symetry of the sign of A is indifferent. So we have found only one A > which is suitable for 
this face. Next we verify that it is also suitable for the others faces. We can then deduce the 20 
vertices of T. 

Now we construct the faces of T. Each of them has five edges which are the shortest geodesic Q 
joining two adjacent vertices Si and Sj. The geodesies can be described as follows [16]. A path I 
on is a geodesic if and only if there is a 2-dimensional plane 11 in passing through the origin 
such that / C n n 5^. All the geodesies are circles. 



(III.l) 



q^G 3a > 0, 3/3 > q = aOSi + f30Sj and \\q\\ = 1 

3a > 0, 3/3 > q = aOSi + (30Sj and 1 = a^ + /J^ + craf3. 



Then a face is the set of shortest geodesies joining two points of the edges. One may also say 
that a face is the set of the projection on of all barycenters of its five vertices in M^. We note 
Fj' the set of all barycenters in of the vertices of Fi. It is included in a 2-plan of equations 
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Q = and ax' + by' + cz' = ^ with even permutations of {a,b,c) = (±-,±1,0). As Fj 

V(xo, X, y, z) G Fj, ax + hy + cz 



we deduce: 

Xq 



(See Appendix A for a detailed description of the faces). The proof is complete. 

Proposition III. 2. The set of t{Si) for all t in I* , and all Si of T has 600 vertices given by: 

(1) a set of 24 vertices given by ^;^(±2, ±2, 0, 0) and all its permutations, 

(2) a set of 64 vertices given by ^^(±\/5, ±1, ±1, ±1) and all its permutations, 

(3) a set of 64 vertices given by ^^(±cr, ±0", ±o", ±^) and all its permutations, 

(4) a set of 64 vertices given by ^^(±o"^, ±^, ±^, ±^) and all its permutations, 

(5) a set of 96 vertices given by ^^(a^, ±^, 0, ±1) and all its even permutations, 

(6) a set of 96 vertices given by ^;^(±\/5, ±^, 0, ±o") and all its even permutations, 

(7) a set of 192 vertices given by ^;^(±2, ±1, ±^, ±0") and all its even permutations. 

They are the vertices of 120 regular dodecahedra which tesselate S^. 

This result is obtained by an explicit calculus of t{Si) for all t in I* and all vertices Si of J^. 
Anyone of these regular dodecahedra is a fundamental domain. Following figure [T] shows the three 
last coordinates of these 600 vertices, viewed from a face of the centered dodecahedron. Straight 
lines between two vertices symbolize edges of pentagonal faces. We note that these coordinates 




Figure 1. 120 CeU 



look like those of Coxeter [6], up to an odd permutation of the three last coordinates. They are 
adapted to I* , unlike those of j6]. From proposition III.l we deduce the following: 
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Proposition III. 3. p{J~), denoted T^, is a centered regular dodecahedron in such that 

(x, y, z) G - 2;2 - y2 _ ^'^^ z) € F. 

It is included in the first ball of the visualization of S^. Its 20 vertices are ^;^(ib^, 0, ±1), 



Each face Fi ^ of F„ is a regular pentagon included in an ellipsoid 



Fi^y C {{x, y, z) e M^, cj^ {ax + by + cz)"^ = 1 - - 

with {ai,bi,Ci) equal to (±^,±1,0), up to an even permutation. 

The set of all barycenters in of the vertices of Fi^y, denoted by Fj'^, is included in a 2-plan of 
equation aiX + hy + CiZ = 

See Appendix A for a detailed description of F^. The following figure is a diagramm, and not a 
visualization, of F (or Fy) because their faces are not in a plan of (or M^). 




Sl4 Sl2 



Si k 




Figure 2. Faces Fi f, viewed from Fi The dashed lines are hidden. 



IV. Equivalence relations on F and F^ 



In order to have relation II. 1 the equivalence relation must identify any pentagonal face of F 
with its opposite face, after rotating by ^ in the clockwise direction around the outgoing axis 
orthogonal to this last face. First specify our notations. We consider the Clifford translations gt 
that have been used for the construction of F, with: 

1 ii 



91 

94 •— 2" ^ 2 o-^ ^ 2J' ya •— 2"^-^ ^ 2 o-J ^ 2"-' yo •— 2 " ^ 2 ^ ^ 2 cr ' 

They are such that: Vi G {1, 6}, gi{Fi) = -Fi+i. The inverses of these six first translations are 
the six other translations: 

_ 1 1 

2" " 



510 := (54)^^ = \al + \li-\i, 



1; , 1 I1, 

o(T1 - ol + 2^k, 



'1 



_ 1 1: _ It, 
oCr± 2 o-J ' 



1^1 1 

2 - 2 

512 := (56) ^ = ocrl 



59 := (53 



58 := (52) 

2"-^ ' 2cr* 2J' 511 (55) —2"-^ 2o-J 2"-' nrz ■— \HK)) 

They all have a translation distance x equal to |. We have for the vertices (see figure [2] for nota- 
tions) : 

5i('S'3) = 'S'e, 5i(5'i8) = 'S's, 51(5*16) = •S'n, 9\{S^) = Svj, gi{S2o) = S2, 
52(5*18) = 5*15, 52(5*12) = 5i7, 52(59) = 52, 52(57) = 5i9, 52(53) = ^4, 
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9'3('S'3) = 5'i, 53(5'7) = 5*11, g3{Sio) = Si7, gsiSu) = Si5, 93{S2o) = S13, 
9i{S2o) = Sq, 54(514) = 5*8, g^iSiQ) = Sii, giiS^) = Si, 94(35) = S12, 

95(35) = Sio, 95(84) = Sq, 95(5*15) = 5*8, 55('S'l3) = 5*9, 55('S'l6) = S7, 
56(5*16) = Sig, 96(Si3) = S2, 56(5*1) = 5*6, 56(5*12) = 5*10, 56(5*18) = 5'i4. 

We define the relation ~ by specifying the equivalence classes q of any q G J^: 

yqeT, q:=I*({q})nT. 

J- has been constructed such that | is the translation distance beetwen two opposite faces. Oth- 
erwise ? is the smallest translation distance of the elements of I*. So: 

5 

q£T^ g = I gi{q), i G {1, 12}} n T. 

It follows that: 

o 

• If q belongs to J^, then q has only one element, 

• If (7 is a vertex of J-, then q has four elements, 

• If q belongs to an edge of a face, without beeing a vertex, then q has three elements, 

• If q belongs to a face and does not belong to an edge, then q has two elements. 

So we have: 

Proposition IV. 1. We define the equivalence classes on T by: 
q(^T ^q = {q}, 

^(hj,k), qGFiflFjnFk ^ q = {q,9i(Q),9j(Q),9k(q)}, 
(^(iJ),<l^Fi<^Fj) and (yk^i,j, q i Fk) =^ q = {q,9i(Q),9j(Q)}, 
(3i, q G Fi) and (Vj ^ q ^ Fj) ^ q = {q,9i(q)}- 

Here the integers k belong to {1, 12}. 

The equivalence relation ~ on F^ is easily deduced from this one on F. Let us denote gi^v the 
application — induced by gi on F^. Hence for all i in {1,...,6}, gi^v(Fi^v) = -^i+i.t; and 
5i> = gi+6,v We have: 

Proposition IV. 2. We define the equivalence classes on Fv by: 

X eFy = {X}, 

3(i,j,k) X £ Fi^^ n Fj,^ n Fk,v ^ X = {X, gi^y(X), gj^^(X), gk,v(X)}, 

(3(i,j), X G Fi,, n Fj>) and (Vk / i, j, X ^ Fk,v) ^X = {X,gi,^,(X),gj,,(X)}, 

(3i, X G Fi^y) and (Vj i, X ^ Fj,^) ^X = {X, g^,v(X)}. 

Here the integers k belong to {1, 12}. 

The geometrical meaning of this equivalence relation is clear : we identify any pentagonal face 
of Fy with its opposite face, after rotating by | in the clockwise direction around the outgoing axis 
orthogonal to this last face (see Figure [2]) . 

V. The Wave Propagation on the Dodecahedral Space. 
We consider the lorentzian manifold Mt x S'^/I* endowed with the metric 

g^j_udx^dx^ = dt^ — ds^^, 
and we study the scalar covariant wave equation associated to this metric : 
(V.l) d"^"^ - Ak^ = 0. 
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Here Ak is the Laplace Beltrami operator on K, that is defined by 

Ak := ^9^5^^ ^5,, g-^ = {g^ , \g\ -- 
V Ifl 



det 5^1,1 



Since K is a smooth compact manifold without boundary, — Ak endowed with its natural domain 
{u E L^(K); Aicu G L^(K)} is a densely defined, positive, self-adjoint operator on L^(K) and the 
global Cauchy problem is well posed by the spectral functional calculus : 

nt) = cos (ty=AK) ^(0) + ''""^'/z^h tm- 



V^Ai^ 

We shall use the functional framework of the finite energy spaces. Given m G N, we introduce the 
Sobolev space 

if™(K) := {u G l2(K), V^u G ^^(K), \ a \< m} . 

where Vk are the covariant derivatives. We can also interpret this space as the set of the distribu- 
tions u G H"^(S^) such that u o g = u for any g ^ I* . Then the standard spectral theory assures 
that for ah '^o G -^^K), ^-i G ^^(K), there exists a unique ^ G (M+; H^(K))nC^ (Mt+; L^(K)) 
solution of (V.l) satisfying 

(V.2) 



^'(t = o) = ^'o, dt^{t = o) = ^i, 



and we have 



/ I dt^it) |2 + I Vk^-W P dfiK = Cst. 
Jk 

Also we have a result of regularity : when G H^iK), G H^(K), then ^' G C° {Rf; H'^iK)) n 

L'(K)). 



H 1 - 



Obviously ^ is entirely determined on x K by its restriction to Mt x J~- To perform the numer- 
ical computation of this solution, we take the domain of visualization C M'^ of the fundamental 
polygon C 5^ as the domain of calculus. Therefore we introduce the map / that is one-to-one 
from Tv onto T defined by 

X = {x,y,z) ^ (v/1 - a;2 - y2 _ z'^,x,y,z) = {fi{x,y, z), f2[x,y, z), h{x,y, z), fi{x,y, z)) , 
and we put 

ip{t,x,y,z) := '^{t,f{x,y,z)). 
We introduce the usual Sobolev space H"^ for the euclidean metric of M'^: 

H"\T^) = {u e L\T^), VaGN^ \a\<m, d" ,u £ (Ty)} . 



Proposition V.l. ^' G C° (M+;F1(K)) n (M+;L2(K)) is solution of (V.l) iff il^ belongs to 

C° [Rt;H^{T^)) n [Rt;L^{Tv)) and satisfies the equation 

(V.3) att^-Aj-> = 0, {t,x,y,z)€R+ xj^^, 

where 

.V4. = {l-x'-y'-z')dl + {l-y')dl + {l-z')dl 

^ ' -2xy - 2xz df^ - 2yz Sfg -3xdi- 3y 82 - 3z 83, 

and the boundary conditions 

(V.5) V(t,X,X') G M X X X ^ X' ^%b{t,X) = ^{t,X'). 

Moreover ^ e (M+;/i'2(K)) n {R+;H^K)) D (M+;L2(K)) iff ^ e (M+; //^(j;)) n 

Ci (M+; H\F,)) n C2 (M+; L^{F,)) . 
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Proof. We denote g the metric induced on J^y by the metric of S"^ and the map /. We 
the coefficients of gij are: 



9 

z,et that 



4 

i=i 



So 



xz 



4 

■E 



1— ?/2_2;2 



1— j/2_ 
1— X'^— J/*^- 



1 — X^— 1/*^ 



1*J 




and 



det (7 



1 



2-2 _ y2 



.2- 



TTT^itS'^'^vlsl^;^ obtain the expression of Ajr 



As Aj-„ : 

Given u G H^{J-v), the trace of n on is weU defined since the domain F.^ is Lipschitz and C° 



piecewise. ij: satisfies the boundary conditions (V.5) iff its pull-back ^' on can be extended in a 



solution defined on the whole Poincare dodecahedron K. The extension to the smooth solutions is 
straightforward. The proof is completed. 

To handle the boundary when applying the finite element method, it is very convenient to take 



into account the boundary condition (V.5 ) by a suitable choice of the functional space. We introduce 



the spaces W^{Fy) that correspond to the spaces H"^(K) : 
1 < m, W^'iTy 



L^{Tv,{l - - - z^) '^dxdydz), 

{u G H'^iJ^y), y{X,X') e dJ^l X r^X' ^ u{X) = u{X')} , 



endowed with the norm 



E 11^" 

\a\<m 



U 



In particular, we have 



W\Ti,) = {ue H^{T,a), X ^ X' ^ u{X) = u{X')} 



and 



The numerical method to solve the Cauchy problem will be based on its variational formulation. 

Theorem V.2. Given ipo £ W^{Tv), ipi £ W^{Tv), there exists a unique 

(R^; W^{J='y)) n C2 (m+; solution of the equation Fv^, and satisfying 



(V.6) 



V^(0,.)=Vo(.), dtm-) = M-)- 
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ip is the unique function in C° (Rt;W^{T^)) n {Rt;W^{J'v)) n (Mt+; satisfying 
>r any (p £ 

{1 — — y^ — z'^)^2'ip(^t, X, y, z)(l){x, y, z) dx dy dz 



df^ 

(V.7) +/ {1 - x'^ - y'^ - z^)-iV'tlj{t,x,y,z) ■V(f>{x,y,z)dxdydz 



V.6) and such that for any (p £ W^{J-y), we have : 
d^ 

2 „ 2 

\yL — X 

{\-x^ -y^ - z^Y^ [(x, y, z) ■ Vip{t, x, y, z)] [(x, y, z) ■ V</>(x, y, z)] dx dy dz = 0. 

Proof: The existence and the uniqueness of the solution of the Cauchy problem are given by the 
previous proposition since the boundary conditions are imposed by our choice of space W^{J-y). Now 
the mixed problem can be expressed as a variational problem. -0 is solution iff for all (j) G W'^{Fy), 
we have : 

/■ 1 

(V.8) =^ /^^ y^=;fryfr^^(*'^'^'^)'^(^'^'^)'^^^2/dz 

{^J^vi^){t,x,y,z)(j){x,y,z)dxdydz. 



\/\ — x'^ — y'^ — z^ 



To get the formulation (V.7), the key point consists in expressing in a symetric manner the last 
integral, without using an integral on the boundary. To invoke the Green formula, we denote 
v{x^ y, z) the unit outgoing normal at (x, y, z) belonging to face We know explicit form of ^j^^, 
the application — t- induced by gi on T^. Thanks to relations written in Appendix A, we 
have: 

5i,„(x,y,2;) = - ( -a -1 ~^ j ( ^ j ' f2.^(^' ^) = 2 ( ~^ ^ ~W \ ^ 
1/ 0"l-\/x\ -./ - crl\/x 



-1 


1 




_\ 

(7 


a 




a 


i) 


—a 


1 










I 

a 


a 




a 


-1 




-1 


_i 






(T 




-1 


_1 


—a 


1 

a 


a 




a 


-1 


—a 


1 





g3,v{x,y,z) = -1 ^ a j \^y j , 54,i,(a;, y, z) = - |^ a -1 ^ j yy 
g5,v{x,y,z) = ^ [ 1 I -a ] I y ] , g6,v{x,y,z) = \i ^ cr -1 ) ( y 

V -1 / V ^ / V 1 ^ / V ^ 

Moreover we know zv(x, y, z) thanks to the equation of the ellipsoids (Appendix A). We verify that: 

Since (n o gi^^p. ^ = u\p. ^, we have for u S W'^{Fv) 

du{X)u{x,y,z) = giAiy{x,y,z)].Vu{giAx,y,z)) = -d^(^g^^^(^^^y^^))u{giAx,y,z)). 
We deduce that for u £ W'^ {Fy) , v £ {F^) , we have 



v{x, y, z)dy(^^^y^^)U da{x, y,z) = - v{x, y, z)d^^^^y^^)U da{x, y, z) 

v{x, y, z)d^i^^^y^^)U da{x, y, z), 
and therefore 

v{x, y, z)di,ix y^)U da{x, y, z) = 0, 
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where da denotes the volume measure on {dTv,g)- We conclude that : 

/ (Aj-^V) 2:, y, z) (p{x, y, z) y/det g dxdy dz = — I g{grad, ijj, grad (f)) 1/ dct g dx dy dz. 

To simplify the writing we shall note in the following X instead of (x,y,z) and |Xp instead of 
+ + for (x, y, z) G F^'-, and dX will designate dx dy dz. Therefore we compute : 

/ (l-|X|2)-lA^Xt,X)0(X)dX = 
- j [(l-x2)9iV(i,^)9i X) + (l-y2)a2V'(i,X)92</'(i,X) + (l-z2)a3V'(i,X)930(i,X)] dX 

+ / [ xy (di^it, X)d2<t>{t, X) + d2^{t, X)di(f){t, X)) + xz {di^{t, X)d3<f>it, X) + ds^it, X)di<t>{t, X)) 

+ yz {d2^j{t, x)d34>it, X) + dsiPit, x)d24>it, x))] dx. 

And finally we calculate : 

/ {l-\X\Y'^A^^^{t,X)4>iX)dX = 

- I ^(t, X) ■ V(/)(t, X)dX+ I (1 - (X • Vi){t, X)) {X ■ V(l){t, X)) dX. 

The proof of the theorem is complete. 

We solve this variational problem by the usual way. We take a family Vh, < h < ho, of finite 
dimensional vector subspaces of W^{Tv)- We assume that 

^o<h<hoVh = W\Tv). 
We choose sequences ^o,/i) V'l,^ ^ such that 

V'0,7i -> ■00 in W^{J'v), V'i,?i -> "01 in L^iJ^v)- 
We consider the solution il^h G C°°(Mt; Vh) of 

V</.ft eVh, ^ (1 - iXfrhhit, X)ct>h{X)dX + ^ (1 - |X|2)-§ V^ft(t, X) ■ VMt, X) dX 

- [ {i-\x\YHx ■'7Mt,x)){x ■vMt,x)) dx = 0, 

satisfying il^h{^, •) = V'o,fe(-)' 9tijjh{0, •) = V'i,/i(-)- Thanks to the conservation of the energy, 

/ (1 dtil^{t,x,y,z) p + I VT^ip{t,x,y,z) |^) ^/ det g dx dy dz = Cst, 



this scheme is stable 



vr>o, sup sxxp Uh{t)\\w^ + \\'^Mt)\\L-^ <Qo. 

0<h<ho 0<t<T at 



Moreover, when tp e C"^ {^t'^ W^^(-^t)))) it is also converging : 

vr>o, sup \\Mt)-Ht)\\w^ + \\^Mi)-TM^^\\L^^^^ ^^0- 

o<t<T at at 

If we take a basis (ej)j^^^.^^ of Vh, we expand iph on this basis : 
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and we introduce 

uit) := *w,V'2^••• .v-^vj 

^=i^ij)l<^,J<N,^ ^=iDij)l<i,j<N,^ J^=(^^.-)l<,J<Ar, 
" v 1 — 1^1 

" v — 1^1 

Jj^v dl — \X\^ 

Then the variational formulation is equivalent to 

MX" + (K + D) X = 0. 
This differential system is solved very simply by iteration by solving 

yi^ljn+l _ 2ljn ^ ^ ^ ^ q 

We know that this scheme is stable, and so convergent by the Lax theorem, when 

< {K + B)U,U > 4 



!l% <MU,U> 
Therefore if there exists K > such that 

V 



< 



AT2' 



the CFL condition 



x,y,zYh 



{l-\X\ 



< 



K 
~h 



(ph 



(1- |X|2)4 



KAT < \/2h, 

is sufficient to assure the stability and the convergence of our scheme. 

VI. Numerical resolution 

VI. 1. Mesh. Our goal is to built a mesh of dT such that it has a fixed size in advance, all the 

edges are splitting in the same way, and the meshes of two opposite faces arc the image from each 
other by the Clifford translation gi fitting for ~. Of course we also want that all points of the mesh 
of Fi^y C dT are on the suitable ellipsoid, and that all points of the edges of Fj^^ are on the right 
geodesic. 

First of all we construct the boundary dJ-y . Each is included in a 2-plan of R'^ and can be easily 
meshed with a convenient metric. We choose to built a mesh of F^, with Fi := (5*3, Sig, Siq, S5, S2o)- 
We denote F^^ := p{Fi) the visualization of F^. As we want the vertices on each edge of Fi are 
equidistant, we choose a priori the desired spherical distance between two consecutive vertices. 
This determines the number of vertices on an edge. Then, by a suitable application, we map F^ in 
the 2-plan 2; = of M'^, which allows us to use a 2-D mesh generator able to respect a given metric 
matrix M. Let us recall the inclusions 
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We do a translation t of Ff^^ by the vector OM5^20 = (0' 0) where M5^20 denotes the middle of 

'S'5S'2o, followed by a rotation r in M'^, with an angle — ^ and axis u = SigM^ 20 = ( at — ; 0)- 

So, thanks to the quaternionic calculus: 



r (xi + yj + -zk) 



1 

V2 



/3-(7 



V2. 



2 crv^3^- 



[xi + yj + zk] + 



'3-cr 



V2_ 



/3-0-' 



and, 



3-cr 



r(x,y,z) 

If (x, y, z) belongs to F^^ we get 

(r o {x,y,z) - 



cr{Z-a-)y + 

1 



/3-(7 



'(7(3-(7)^ + tT2(3-a)y + 



/3-cr 



/3-(T' 



a 



V3 



a 



y + 



2V2 



(x',y',0). 



Therefore r o t{Fi^) is in the 2-plan z = of Of course we want a metric on M endowed by 
the metric of 5'^. We denote by / the application that maps r o t{Ff^) to Fi, that is / : 



rot{Fl,) 

{x,y) 



1-^ 



{x,y,0) 



followed by: 



F^ C 



Ft ^ 
(x',y',z') := (rot)-i(x,y,0) ^ 



F^ C 



,x ,y ,z 



2V2 



x',y' ,z' 1 1 



,x ,y ,z 



=1,4 



We deduce that the coefficients of M are given by: 
(VI.l) ma ■■= - fi y) ^iifj ' : 

See Appendix B for the values rriij. 

We used the software FreeFem++ ([8|) for the generation of the 2-D mesh of r o t [Fi^). Then, 
applying /, we get a mesh of Fi with the wanted size. Moreover, as we know the equation of the 
2-plan containing Fi^, we are able to control and to correct the coordinates of the vertices of the 
mesh of F^^ with a great precision. We can also control and correct that the vertices on the edges 
of Fi are on the convenient geodesic (III.l), distant from each other with the wanted distance. 

Then, by the use of rotations in M"^ we get a mesh of the five adjacent faces FI j,, , Fg ^ (see 

Appendix C). After normalizing we have got a mesh of Fi, Fq. At last, by the use of gi, 
qq we deduce a mesh of F7, F12. Once more we control and correct the vertices on edges. As 
we also know the equation of the ellipsoid that contains each face, we also control all the vertices. 
Meshes of each Fj^„ are obtained from those of Fj by discarding the first coordinate of the vertices. 

o 

To get the mesh of F we used the software Tetgen ([18j) in several ways. We have created a 
mesh with an inner sphere or an inner dodecahedron, in order to impose two different average sizes 
of tetrahedra depending on the location. We also used Tetgen with a .mtr file that gives us the 
wanted size near each point. In all cases we employed an option so that there is no point introduced 
into dTv 
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We can see the accuracy of our mesh by calculating the sum of the volume of tetrahedra which 
compose This sum must be closed to ff^ that is ~ 0.1644934066. 



Mesh 


number of vertices 


number of nodes 


number of tetrahedra 


sum 


relative error 


F089 
FlOl 


182162 
233858 


730309 
1600118 


4642744 
10198710 


0.1644927119 
0.1644928648 


4,22^-006 
3,29^-006 



The minimal and maximal distances beetween two vertices of the bounder of these meshes are: 



Mesh 


minimal 


maximal 


F089 
FlOl 


2.06^ - 003 
1.82^; - 003 


6.11^;- 003 
4.90£; - 003 



We can also look at the tilling of 5^ . In figure [s] we display six images of Tv , as well as J- in the 
first ball of visualization of . All points of theses images have a positive first coordinate in M^, 
except those of g2 gi gi{T)- To fully represent this last set, we should also draw the second ball of 
visualization of 5^, in which its representation is the same as in the first ball. 




MFv, Mgi{J^v), Mgi gi{J^v), ■56(-^d), Mgegei^y), ""g^ geiJ^v), ■5251 gii^v) 
Figure 3. The Fundamental Domain (dark blue) and six images of it. 



VI. 2. Vh space. We cons truct the finite element spaces Vh of Pi type. We take into account the 
boundary condition (V.5) in the definition of the finite elements, so that Vh C W^{J-y). We note 
Th all tetraedra of a mesh, and h 



Vh := {v : 



UxeTh^- Then we introduce : 
V G C°(j;,;,), VK G Th, v\K G IPiW, M^M' ^ v{M) = v{M')] . 



The equivalent points on dJ-v are known thanks our construction of meshes. The number of nodes 
Nh is the sum of x ^ nbve + ^ x nbvf + nbvi + with nbve the number of vertices on an 
edge of a face that are not a vertex of Tv, nbvf the number of vertices on a face that are not on 

o 

an edge, nbvi the number of vertices in T^. 

If j is the number of a node and if Mj denotes a vertex of the mesh, we construct a basis 
of Vh by: 



(1) If J is associated to a node that does not belong to dTy : e^{Mi 
There are nbvi functions of this kind. 
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(2) If j is associated to a node that is a vertex of J-'^: 

■' \^ U otherwise. 

There are five functions of this kind. 

(3) If j is associated to a node that belongs to a face of dJ-^ and not to an edge: 
hr.^A_ \ 1 if Mi = Pi, 



1 otherwise. 
There are 6 x nhvf functions of this kind. 
(4) If j is associated to a node that belongs to an edge of a face of dT^ and is not a vertex of 

^ ' [0 otherwise. 
There are 10 x nhve functions of this last kind. 



VI. 3. Matrix form of the problem. K.(z,j), D(z,j) and M(i,j) are found with a numerical 
integration using CUBPACK These matrices are sparse and symetric. So we choose a Morse 

storage of their lower part, and all of the calculations will be performed with this storage. To solve 
the linear problem we use a preconditioned conjugate gradient method. The preconditioner is an 
incomplete Choleski factorisation, and the starting point is the solution obtained with a diagonal 
preconditioner. 



VI. 4. Initial data. We choose different initial data. Some of them have a more or less small 
support near given points, others are chosen randomly. For the wave depicted in figure |4j we have 
taken 

V'o(^) = lOOCf-^'-'^ol'-'-o, for d(X,Xo) < ro, and V'o(^) = 0, for d(X,Xo) > tq. 

with = (0, 0, 0) and ro = 0.3. In the following we refer to as the name Initc- And for the 
one depicted in figure [6j we have taken a similar function with a smaller support, and especially a 
support not centered at the origin. In the following we refer to as the name Initexc- 



VI. 5. Time resolution. In order to see the stability of our method, we compute E^it), the dis- 
cretized energy at the time t: 

E{nAt) = (M — , \ + ((K + D) X'-^ , X"> . 

It is well known that our scheme is conservative, hence Eg^ must be invariant all along the resolution: 
this is the case. For example, with the previous initial data and the use of the mesh F089, we obtain: 



Initial data 


EM 


Edm 


Initc 
Initcxc 


9328.34990491949 
6248.87928527422 


9328.34990491958 
6248.87928527418 



Here are some characteristic pictures for a centered initial data. At each time there are a view of 
the solution on dT^j, and on a cutting plane {z = 0). 
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« 




□ 





Figure 4. A transient wave at t = 0, t = 2.8, t = 5.44, t = 11.52, t = 57.92, t = 74.24. 

The solution varies extremely fast as you can see on the next picture. This was expected since 
eignenvalues are large. 




Figure 5. The solution -'^o) from t = to t = 80. 




SI « 



a B H H H H 



Figure 6. Another transient wave at t = 0, t = 2.65, t = 14.65, t = 24.5, t = 33.45, 
t = 71.10. 
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Some movies are available on this personal home-page 



VI. 6. Eigenvalues. We test our scheme in the time domain, by looking for the eigenvalues that 
are explicitly known (see [T], [9], [10], [H]). Since the Laplace-Beltrami operator Ak on PDS is 
a non positive, self-adjoint elliptic operator on a compact manifold, its spectrum is a discrete set 
of eigenvalues —q^ < 0, and by the Hilbert-Schmidt theorem, there exists an orthonormal basis in 
L^(K), formed of eigenfunctions {'>pg)q C i/°°(K) associated to g^, i.e. 

One has: 

with /? G {1, 13, 21, 25, 31, 33, 37, 41, 43, 45, 49, 51, 53, 55, 57} U {2n + 1, n > 30}, and we take tpo = 
Therefore any finite energy solution ^|;{t, X) of dfi/j — Aj^ip = has an expansion of the form 

^qe'^'^^ipq{X). More precisely, if we denote <, > the scalar product in L^(K), we write 

ij{t,x) =^(<atV(o,.),i >t+ < v(o,.),i >) 



To compute the eigenvalues q'^ we investigate the Fourier transform in time of the signal il){t,Xo) 
in the case where dtip{0,X) = 0, for different Xq. Practically, during the time resolution of the 
equation, we store the values of the solution at some vertices M{X) for the discrete time kAt, Ni < 
k < Nf. We choose the initial step Ni in order to the transient wave is stabilized, that to say NiAt 
is at least greater than the diameter of J^y, i.e. NiAt > 2 * d{0, Si) = d{Si, Su) = arccos ( ^"^g~^ ) — 
0.776279. Then we compute a DFT of {iph{kAt, ^))Arj<fc<Arj, with the free FFT library fftw3. After 
we search the values jmaxi, jmax2, ■■■ for which the previous result (^j(X))o,Arj^-Arj+i has a 
maximum. The eigenvalues found by the algorithm are expressed as: 

q = ——imax. 

^ {Nf -Ni + l)Ar 

We have made tests by varying parameters such as: Ni, mesh, initial data. In almost all cases we 
have obtained the 36 first consecutive eigenvalues with an relative error of about 10~^. For example 
here are non zero values found with the mesh F089, Initexc as initial data, and At = 0.0005. FFT 
is performed from iV, = 1800 to Nj = 106770. 



http: //www. math. u-bordeauxl . f r/~mbachelo/PDS .html 
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/3 


exact q 


numerical result 


relative error 




/3 


exact q 


numerical result 


relative error 




13 


168 


167.6126 


2.3060753E' - 


03 


69 


4760 


4871.141 


2.3348872E' - 


02 


21 


440 


438.8956 


2.5100708E' - 


03 


71 


5040 


5171.347 


2.6060946i!J - 


02 


25 


624 


628.4641 


7.1539269£/ - 


03 


73 


5328 


5463.117 


2.5359832£^ — 


02 


31 


960 


965.4485 


5.6755068ii^ - 


03 


75 


5624 


5762.894 


2.4696665E' - 


02 


33 


1088 


1093.790 


5.3216149£^ - 


03 


77 


5928 


6089.032 


2.7164599£; - 


02 


37 


1368 


1374.492 


4.7458406£^ — 


03 


79 


6240 


6424.146 


2.9510655E - 


02 


41 


1680 


1696.904 


1.0062009E' - 


02 


81 


6560 


6787.617 


3.4697667E' - 


02 


43 


1848 


1865.750 


9.6047800E' - 


03 


83 


6888 


7121.304 


3.3871040E' - 


02 


45 


2024 


2042.602 


9.1906162E' - 


03 


85 


7224 


7462.999 


3.3083960£; - 


02 


49 


2400 


2431.920 


1.3299967£; - 


02 


87 


7568 


7833.520 


3.5084505£^ - 


02 


51 


2600 


2633.263 


1.2793344i5; - 


02 


89 


7920 


8213.017 


3.6997046E' - 


02 


53 


2808 


2842.612 


1.2326320£; - 


02 


91 


8280 


8536.121 


3.0932499E - 


02 


00 




ouoy.yoy 




UZ 


yo 


oD4to 


SRfll /IQfl 

oDUi.^iyu 


o.oiouyou-C/ — 


uo 


57 


3248 


3298.837 


1.5651835^ - 


02 


95 


9024 


8645.207 


4.1976172^;- 


02 


61 


3720 


3789.000 


1.8548321^ - 


02 


97 


9408 


8954.335 


4.8221201-E;- 


02 


63 


3968 


4039.324 


1.7974915^; - 


02 


99 


9800 


9021.282 


7.9461001^; - 


02 


65 


4224 


4297.656 


1.7437559S - 


02 


101 


10200 


9405.365 


7.7905372E - 


02 


67 


4488 


4579.910 


2.0479091E: - 


02 


103 


10608 


9797.453 


7.6409020£; - 


02 



VI. 7. Conclusion. Wc note a good agreement of the eigenvalues obtained by the spectral analysis 
of the transient waves that we have computed. We emphasize that this result is a strong evidence 
of the accuracy of our numerical approximation since the waves are rapidly oscillating due to the 
large value of the eigenvalues. Therefore we conclude that we have validated this computational 
method of the waves on the Poincare dodecahedral space, and we plan to extend it, in the future, 
to the case of the non-linear waves on this manifold. 
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19 



(1) Coordinates of the vertices Si of T: 



52 = 




^2 


5*6 = 


" 171 ( 




5*10 


_ 1 

2^f2 


(- 


5*14 


_ 1 
2s/2 


(- 


'S'ls 


_ 1 


(- 



^9 = ^ ("';-^;^' ^) ' = 171 f '0'/); = 171 1' -^) ' = ^ -1, ^,0) ,^ 

= ^ (f^^-^,0,-l) , 5i4 = ^(cj2,^,-^,^), 5i5 = ^ ((T^^,0,-l) , 5i6 = ^ ((72,-^,-^,-^ , 
^i^ = 17l(^''^'^'-^)' = ^(^^ -1,-^,0), 5i9 = ^(aM,-^,0), 52o = ^(a^0,-l,J.). 



(2) Images by the Clifford translation gi of the face Fj of and of its edges. 
g\ maps F\ to F7, and we have for the vertices and the edges: 

ffi(5'3j_= 'S'e,^ 9]{S\%) = S^^ giiSie)^ = Su, gijSr,) = Sij, gjjSw) = 5*2, ^ 

5'i(5'3'S'i8) = SqSs, giiSigSie) = SsSu, giiSieS^) = SuSn, gi{S5S2o) = S17S2, 
gi{S2oSs) = S2SQ. 

g2 maps F2 to Fg, and we have for the vertices and the edges: 

g2{ Si8) = Si5,^ Jj2iSl2) = 517, 52(59) = S2j _ 52(57) =519, jMSs) =J±^ 
(72(5i85i2) = 5i55i7, 52(51752) = 585ii, 52(5957) = 525i9, 52(5753) = 5i954, 
52(53518) = 545i5. 

g'3 maps F3 to F9, and we have for the vertices and the edges: 

53(53)_=5i,^ 53(57) =_5ii, 53^(5io) = 5i7, 53(5i4^ = 5i5, gjiS2o) = 813, 

53(5357) = 5i5ii, 53(57510) = 5ii5i7, 53(5io5i4) = 5i75i5, 53(514520) = 5i55i3, 
53(52053) = 5i35i. 

g4 maps F4 to Fio, and we have for the vertices and the edges: 
54(520) = 59, (74(514) = 58, 54(519) = 5ii, 54(54) = 5i, 54(55) = 5i2, 
54(52o5i4) = 5958, g'4(5i45i9) = 585ii, 54(51954) = 5ii5i, (74(5455) = 5i5i2, 
54(55520) = SuSg. 

55 maps F5 to Fii, and we have for the vertices and the edges: 

55(55) = 5io, 55(54) = 56, 55(515) = 58, 55(513) = 59, 55(5 *16) = 57, 
55(5554) = 5io56, 55(54515) = 5658, 55(5i55i3) = 5859, 55(5i35i6) = 5957, 
55(51655) = 575io 

56 niaps ^6 to F12, and we have for the vertices and the edges: 

56( 5i6) = 5i9^ _56(5i3)^_52, 56(5i) = 56, _ jf6(5i2)_=^5io, ^jieiSis) =_5i4^ 
56(5i65i3) = 5i952 , 56(5i35i) = 5256, 56(5i5i2) = 565io, 56(5i25i8) = 5io5i4, 

56(5l85i6) = 5i45i9. 
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(3) F-', the set of all barycenters in of the vertices of Fi, is included in a 2-plan of M''. There 
are the equations of these 2-plans. 



Ft 






X, 


y, 




e I 




Xq 


Ft 




(2;o 


X, 


y, 


^) 


e I 
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(4) So, after having normalized the points of F^ ^ we get that Fi is included in an hyperplan of 



There are the equations of these 3-plans. 
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(5) We deduce from the previous items that F^^, the set of all barycenters in of the vertices 
of F ^, is included in a plan of M"^. There are the equations of these plans. 
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(6) As = 1 — x^ — — z^, we obtain from item 4 that the faces Fi^y of Tv are included in an 
ellipsoid. 

and F-j^y are included in the same ellipsoid {(x, y, z) G M^, (cr + 2)x^ + Sa^y^ + 2;^ + la^xy - 
F2 y and Fg are included in the same ellipsoid {{x,y, z) G M^, Scj^x^ + + (cj + 2)z'^ - 2a^xz = 
F^^v and Fg^y are included in the same ellipsoid {(x, y, z) G M^, x^ + (fi + 2)y^ + Scj^z^ — 2a^yz = 
Fi^y and FiQ^y are included in the same ellipsoid {(x, y, z) G M^, {cr + 2)x^ + Sa'^y'^ + — 2a'^xy 
F^^y and Fn^y are included in the same ellipsoid {(x, y, z) G M'^, x^ + (cr + 2)y^ + Sa^z^ + 2a^yz ■ 
F^^y and Fi2,t; are included in the same ellipsoid {(x, y, z) G M'^, 3cj^x^ + y^ + (ct + 2)z^ + 2a^xz 
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VII. 2. Appendix B. W e give the expression of the metric matrix on endowed by the metric 
of 5^ which is written in VI. 1 We denote by / the application that transfoms the 2-D mesh in a 
mesh of Fi , that is / : 
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We simplify: 




9{x,y) 
202 



Thanks to Maple we obtain: 



X [I6OO x2y2 + 1600 + 80 V^V^x^y + 400 V^x'^y + 160 V^V^xy"^ + 800 V^y^ + 160 V^Vby^ 
+860 x^ + 340 V5x^ + 80xy + 80 \/5xy + 760 V5y^ + 2000 + 86 ^/2^/5x 
+170 \/2x + 258 V2V5y + 610 V^y + 845 + 374 a/s] , 

77122 := 160 [80x2 + 80y2 + 8V2V5a; + 20V2y + 4V2V5y + 45 + 17V5]~^ 

X [32OOX'' + 3200x2y2 + 640^2^5x3 + 800^2x2?/ + 160V2^x2y + 320^/5^/2xy'^ + 3800x2 
+1320v/5x2 + 160V5xy + 160xy + 1680y2 + GAOVEy'^ + 660V2x + 348V2V5x + 244^2^5^ 
+580v^y + 717^ + 1625] , 

77712 := 32000 [80x2 ^ gQy2 ^ 8^V5x + 20^2?/ + 4^V5y + 45 + UVE] ~^ 

X [8192000x^7/ + 32768000x'^7/3 + 49152000x^7/^ + 32768000x3/ + 8192000x7/^ + 204800\/2\/5x*^ 

+1024000^/2x9 + 3686400^/2^/5x8^ + 2457600^/5^/2x^y2 _^ i2288000^/2x'^7/2 ^ 11468800^2^5x^7/3 

+6144000V5\/2x5/ + 30720000^/2x5/ + 12288000V2\/5x^/ + 5734400\/5\/2x3/ 

+28672000V2x3y6 + 4915200\/5\/2x2y'^ + 9216000v^xy8 + 1843200^5^2x^8 + 409600\/5\/2y9 

+921600x8 + 921600^5x8 + 7782400x/5x'^7/ + 27443200x'^?/ + 8601600x^^2 + 8601600^5x^^2 

+84787200x5/ + 25804800 \/5x5/ + 15360000\/5x''7/ + 15360000x^/ + 87244800x3/ 

+28262400^/5x3/ + 8601600x2/ + 8601600^/5x2/ + 2990_0800xy'^ + 10240000^/5x7/^ + 921600^5/ 

+921600/ + 1495040^2 VSx'^ + 3993600x/2x7 + 7884800a/2 1/5x^7/ + 13619200\/2x^7/ 
+13578240^2^5x^7/2 + 35328000^2x^^2 + 1832960072^5x^^3 ^ 32256000\/2x4y3 

+22835200^/2^/5x3/ + 57856000^/2x3/ _^ ii857920\/2/5x2y5 + 21196800\/2x2/ + 26521600\/2xy' 
+10752000\/2^/5xy6 + 1413120^/5^/2^^^ + 2560000^/2^^ + 4802560x6 + 2365440V^x6 + 43054080x5y 
+18201600\/5x5y + 31795200x'^/ + 15513600V5x'^/ + 95078400x3/ + 40704000x/5x3/ 
+15820800\/5x2/ + 32716800x2/ + 51655680x/ + 22379520^5x7/^ + 5232640/ + 2508800\/5/ 
+2740224V^^/5x5 + 6259200^/2x5 + 19347200/2x^7/ + 8878336/2/5x4y + 16392192/2/5x3/ 
+37248000/2x3/ + 27302400/2x27/ + 12486144v^/5x27/ + 13570048/2/5x/ + 30681600/2x/ 
+5241600V^7/5 2389248/2/5/ + 5671424x'^ + 2559744/5x'' + 32037888x3y + 14255616/5x3y 
+9666048/5x2y2 + 21451776x27/2 + 36815872x/ + 16409088/5x/ + 6244864/ + 2809600/5/ 
+1914304/2/5x3 + 4288192/2x3 + 8825088/2x2?/ + 3954816/2/5x2y + 5696448/5/2x7/ 
+12752064/2xy2 + 3657984/2/ + 1638528/2/5/ + 2327712x2 + 1041696/5x2 + 3759032/5xy 
+8408152xy + 1119072/5/ + 2501088/ + 458214/2/5x + 1024706/2x + 357278/2/57/ 
+798798/27/ + 192091 + 85909/5] . 
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VII. 3. Appendix C. We present the rotations that map to each F^^^ for i belonging to {2, 6}. 
Using them, we can built a mesh of the five adjacent faces F2^^, , Fg „ of 



*F^.^ is sent to Fq .^ and F^.^, by rotations in with an angle ib^ and an axis u :- 



where G2 denotes the center of We have: u 
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•Ff^ is sent to ^, and F4 ^ by rotations in M"* with an angle and an axis u : 
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where G3 denotes the center of Fg^. We have: « 



with an angle ^ and an axis u: 
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»Ff is sent to F5 j, by a rotation in M"^ with an angle ±^ and an axis u : 
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Hence : F\.^ — > F^,^. We have: 
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